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Abstract 



In hep-lat/0701018 , Creutz claims that the rooting trick used in simulations of staggered fermions 
to reduce the number of tastes misses key physics whenever the desired theory has an odd number 
of continuum flavors, and uses this argument to call into question the rooting trick in general. Here 
we show that his argument fails as the continuum limit is approached, and therefore does not imply 
any problem for staggered simulations. We also show that the cancellations necessary to restore 
unitarity in physical correlators in the continuum limit are a straightforward consequence of the 
restoration of taste symmetry. 
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1. Simulations of QCD using staggered fermions have, in the last few years, attained im- 
pressive precision in calculations of the hadron spectrum and matrix elements [l|, HI, H, 0, 0] • 
A crucial part of the methodology is the need to use the so-called "fourth-root trick." Each 
flavor of staggered fermion on the lattice represents, in the continuum limit, four degenerate 
continuum fermions (conventionally called "tastes"). To reduce to the desired number of 
continuum fermions — one per flavor — the simulations are done with the fourth root of the 
(four-taste) staggered fermion determinant. Since the staggered Dirac operator does not 
represent four degenerate fermions for a 7^ (a is the lattice spacing), the validity of this 
procedure is a non-trivial issue. In the past year, there has been signiflcant progress on this 
issue, showing both the nature of the discretization effects at non-zero lattice spacing and 
how, given certain technical assumptions, the correct physics is obtained in the continuum 
limit p, 0, H, H, The required assumptions are plausible (at least to us), but have not 



yet been proven analytically or tested in detail numerically. 

This note concerns particular criticisms of the rooting procedure given by Creutz He 
argues that the rooting procedure "mutilates the expected continuum holomorphic behavior 
in the quark masses" and is manifestly flawed when applied to a theory with an odd number 
of flavors. Given the phenomenological signiflcance of the results obtained with staggered 
fermions, it is important to determine the correctness and applicability of his arguments. 
Here we show that his arguments do not apply to rooted staggered fermions as they are 
used in simulations of QCD. In other words, Creutz's criticisms do not call into question the 
correctness of the physical results obtained in the continuum limit of staggered QCD with 
the rooting procedure. 

Before going into details, we stress several points: 

A. Simple, well-tested properties of the staggered theory are all that are required to refute 
Ref. 

B. Showing that Creutz's reasoning is incorrect or misleading does not of course prove 
the correctness of the rooting procedure, but only that his particular criticisms do not 
apply. ^ 

C. We do not dispute that the structure of the rooted staggered theory is complicated 
"ugly" [9I) at non-zero lattice spacing, a. As we have emphasized jsl, Q, S, 0, 10|, there 



are certainly violations of unitarity and locality for a 7^ 0. However, all that is required 
for existing and planned staggered simulations to compute QCD quantities correctly 
is that the rooted staggered theory reproduces QCD in the continuum limit for flxed, 
positive quark masses. Any "mutilations" or problems listed by Creutz at non-zero a 
that go away in the continuum limit are not related to the issue of correctness. Indeed, 
such discretization effects would be relevant at all only to the extent they were not 
understood and therefore not properly taken into account in flnal systematic error 
estimates after extrapolation to the continuum. As explained below, that is not the 
case. 

D. Several properties Creutz flnds "peculiar," "contrived," or "absurd," such as the lack 
of commutativity of an ultraviolet and an infrared limit, or the existence of unitarity 



11 1 does not respond to the arguments of Refs. supporting the 



We note that Creutz's work 
vahdity of rooting, and indeed does not refer to several of these references. 
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violations outside of the physical subspace, are not special to the rooted staggered 
regulator, but can appear with other regulators or in manifestly correct continuum 
versions of QCD. We give examples below. 



2. Creutz's central argument can be summarized as follows. QCD with an odd number 
of degenerate quarks is not invariant under a change in the sign of the (common) quark 
mass, m 12(]. The corresponding rooted staggered theory, by contrast, is invariant under 
this change at any lattice spacing (since it descends from a theory with four tastes per 
flavor, which is itself invariant). Thus, taking the continuum limit, the lattice theory cannot 
correctly represent the desired underlying theory for any m. In particular, odd powers of m 
are missing in the chiral expansion. Creutz claims that this argument is particularly sharp 
at finite volume, because one then avoids the discontinuities that can develop in infinite 
volume. In finite volume, one expects in the continuum limit that physical quantities for 
theories with an odd number of quarks will be analytic functions of m having both even 
and odd powers. Such functions cannot, Creutz argues, be obtained from the rooted lattice 
theory whose partition function is expandable in (only) even powers of m at any lattice 
spacing. 

Two key points are needed to understand the flaws in this argument. The first is that 
rooted staggered fermions correspond, in the continuum limit, to a theory with a positive 
quark mass, irrespective of the sign of the mass in the underlying four-taste quark determi- 
nant. Schematically, 

[det(Atag + m)ff/' > [det(Deont + |m|)]^^ , (1) 

a— >0 

where -Dcont is a continuum Dirac operator for one quark field, and Nf is the number of 
flavors. This equation is schematic as we are not being careful with overall factors or mass 
renormalization, neither of which are relevant here as they do not change the symmetry 
properties under a change in sign of m. However, for those who might worry that problems 
could lurk in this schematic form, we note that it is possible to define precisely the form 
in the continuum limit using renormalization group techniques — see Ref. 0]. Note that 
the root on the left-hand side of Eq. ([1]) is always chosen to be positive.^ Equation ([1]) 
implies that the physical theory with negative quark masses is not directly accessible with 
the rooting trick. 

The second basic point is that rooted staggered fermion determinants are, by construction 
(because of the root), non-analytic functions of the quark mass, even in finite volume. In 
particular, in the continuum limit, they are non-analytic at m = 0. 

We elaborate on these points below, but first show how they refute Creutz's argument. 
This can be seen from Eq. ([1]). In a rooted staggered simulation, the left-hand side (which 



^ Creutz states that "attempts to define the rooting procedure by selecting one fourth of the eigenvectors 
will necessarily involve ambiguities." We agree with this statement, but stress that this is noi what is done 
in the rooting trick, and thus is not relevant to the issue of its correctness. Rooting is an unambiguous 
procedure that takes the unique positive fourth root of det(Z3stag + which is strictly positive for m real 
and non-zero. The effect of rooting is to replace each approximate quartet of eigenvalues by its geometric 
average. Note that when a sufficiently large chemical potential, /i, is turned on, one must deal with the 
(now complex) individual eigenvalues of -Dstag, and as a result ambiguities do arise when taking the fourth 
root of the determinant 



13[. However, this does not apply in the case of = considered here. 
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is the fermion contribution to the weight for a given gauge configuration) is indeed even in 
m, as in Creutz's argument. But the proposed continuum hmit, given by the right-hand 
side, is also even. Thus there is no contradiction. There is no problem with the right-hand 
side being non-analytic in m because the left-hand side is non-analytic by construction. In 
particular, in finite volume, physical observables in the rooted staggered theory can have 
both even and odd powers of \m\ in the continuum limit. 

This argument shows that the rooting procedure can lead to the desired continuum de- 
pendence on m (for m > 0) when a —>■ 0. There is no inconsistency. Whether it does lead to 
the desired dependence is a different issue. This requires that taste symmetry be restored 
in the continuum limit, for which there is, by now, extensive evidence. For example. Fig. 
1 of Ref. [13] shows not only that the splitting of various-taste pions goes to zero as the 
continuum limit is approached, but also that the rate of approach to zero (as Q;|(a)a^) is in 
agreement with the theoretical prediction for the version of the staggered action used. The 
restoration of taste symmetry is also found directly in numerical studies of the eigenvalues, 
which fall into approximate quartets for eigenvalues well below the cutoff, with the splittings 
decreasing toward zero as a is reduced [H, 1^. Further, these studies show that, in topo- 



logically non-trivial configurations, what would be zero modes in the continuum (or with an 
overlap lattice operator) appear as quartets of near-zero staggered eigenvalues. Note that 
there are also good theoretical reasons for believing that taste symmetry is restored in the 



continuum limit. In particular, it has long been known [17|, ll8|| that taste violations are due 
to an operator with dimension greater than four.^ 

We now illustrate how the expected behavior in the continuum limit is approached, 
assuming the restoration of taste symmetry. As stressed by Creutz, the dependence on the 
sign of m in the continuum arises from zero modes of -Dcont- Consider a gauge field with 
winding number z/ = 1, so that the continuum determinant can be formally written as 

[det(D,ont + m)f^' = [mF,,^,{m')f'' , (2) 

where the factor of m comes from the zero mode of -Dcont, and -Fcont('^^) is the contribution 
of the non-zero eigenvalues. Because the latter come in imaginary pairs ±iX, -Fcont(^^) is 
necessarily a function of and is strictly positive for m ^ 0. We are imagining working 
in finite volume so that the spectrum is discrete. For odd Nj the contribution of such 
configurations is manifestly odd in m. 

How does the rooted staggered theory approach this limit? Like the non-zero eigenvalues 
of -Dcont, the eigenvalues of -Dstag are imaginary and occur in complex conjugate pairs. The 
low-lying eigenvalues form into quartets, and these quartets become degenerate in the con- 
tinuum limit. For a configuration with winding number unity there thus will be two pairs 
of eigenmodes with eigenvalues of 0(a), corresponding to one zero mode per taste in the 
continuum limit. Labeling these as ±iaXi^2^\ the rooted staggered determinant is 

[det(Atag + m)ff/' = [{{aXf'^^'^y + m'}{{aXt''')' + m'} F^,.,{m')] ""'^^ . (3) 

Like Fcont(^^), -Pstag('^^) is strictly positive for m 7^ 0. As claimed by Creutz, Eq. ([3]) is 
manifestly an even function of m. If we stay in a finite volume, the determinant is analytic 



The argument has been made more rigorous in Refs. [7|, |9|. Recent calculations have also shown that the 



size of taste violations can be understood quite well perturbatively [19 1. 
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in the complex m plane aside from cuts along the imaginary-m axis. In particular, its 
expansion about m = contains only even powers of m, for a > 0. If we take the continuum 
limit, ^ however, the nearest cuts (due to the near-zero eigenvalues ±iaXi^2^^) move to the 
origin, and, focusing on the zero-mode sector, one finds 

[det(Dstag + m)ff/' [F,tag(m2)] "^^/^ . \mff . (4) 

a— >0 

The key point is that the non-analyticity intrinsic to a root leads to a positive contribution 
from the zero modes, irrespective of the sign of m. In particular, we see how, for odd Nf, 
odd powers of |m| arise in the continuum limit due to this non-analyticity. 

That the rooting trick leads to a continuum theory with a positive mass is not related to 
the use of staggered fermions. Consider the lattice theory obtained by taking the positive 
Nf/4: power of the determinant for four degenerate "tastes" of overlap fermions. There are 
exact zero modes in this theory, and the contribution of each to the determinant is 

= \m\^f . (5) 

The remaining eigenmodes lead to an even function of m. For positive m, the rooting has 
no impact, and one obtains the correct continuum limit, with, for odd Nf, both even and 
odd powers of m. For negative m, one obtains the continuum theory with positive mass |m|. 
All this holds in any finite volume. Thus the rooted overlap theory provides a clear example 
of the failure of Creutz's argument. 

Comparing the rooted staggered and overlap theories, one sees that to reproduce the form 
of the overlap result, it is essential to take the continuum (a — > 0) limit before the chiral 
(m 0) limit. This fact is already well known in the literature [HI, 20, 2l[. Indeed, with 



or without the fourth root, staggered fermions do not always produce correct continuum 
physics if one takes the m ^ limit before the a — limit. ^ Thus the problems that 
Creutz finds in the staggered theory when a quark mass is set to zero before taking the 
continuum limit are irrelevant to simulations, for example those of the MILC collaboration 
dl, S, i, i, 9, in which the quark masses are kept positive for all lattice spacings, and any 
extrapolations are performed in the proper order. 

We note that the effective low energy theory for (rooted or unrooted) staggered quarks, 
staggered chiral perturbation theory |22|, provides an explicit framework to see both that 
the correct continuum physics is obtained, and that the continuum theory is always that 
with positive mass. The chiral theory for rooted staggered quarks, including discretization 
effects, is obtained using a limiting procedure with an even number of flavors x tastes at 
each step. This necessarily results in a theory which is invariant under a change in sign of 
m: for m < one can perform a non-singlet axial rotation to make m positive since there is 
an even number of flavors x tastes. Nevertheless, the effective theory in the rooted case goes 
over as a —i> to the correct, positive quark mass, continuum chiral theory, for both odd 
and even numbers of flavors.^ For example, the masses of pions at leading order take the 
form = 2BQ\m\ + O(a^), with Bq the continuum low-energy constant. The continuum 



^ Which is to be done at_fixed physical volume. 

^ Note that, while Ref. 21 1 argues that the continuum and chiral limits commute for many common ob- 



servables, other observables for which these limits do not commute are also discussed. 
^ At the moment, the effective theory calculations can only be traced from negative to positive masses 
in infinite volume. In finite volume, for small enough masses one necessarily enters the epsilon-regime 
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theory with negative mass (which, for Nf = 3, displays spontaneous CP violation is 
not directly accessible in the chiral theory for rooted staggered quarks, just as at the quark 
level J 

3. Another feature of the rooted staggered theory that bothers Creutz is the existence of 
multiple versions of the physical pions, with various tastes. Indeed, due to the exact non- 
singlet chiral symmetry of staggered fermions in the massless limit, one pion per flavor is a 
true (pseudo-)Goldstone boson, even in situations where no Goldstone particle should exist: 
a one-flavor theory, for example. While Creutz recognizes that this feature may be repaired 
in the continuum limit by cancellations from pions of other tastes, he finds the suggested 
mechanism of repair an "unproven conjecture" that is "rather contrived." On the contrary, 
the mechanism depends only on the restoration of taste symmetry in the continuum, for 
which there is a great deal of evidence, as discussed above. To see how taste symmetry fixes 
this, and many similar "problems," it is necessary simply to examine the continuum-limit 
theory.^ Consider the Euclidean QCD partition function in a formal continuum setup, 

N 

_ I A T^^l. T^ZT. „, — / C I 

-'cont 



/ J 
VAVijVi^ exp(^-Sg + J2i^iiD + mi)ij,^ (6) 
1=1 

= VA exp{-Sg) Y[det{D + mi). (7) 

i=i 

Here Sg{A) is the gauge action, D is the massless Dirac operator, and rrii is the mass of 
the i^^ flavor. We assume all are positive, which implies that the (formal) continuum 
determinant, det{D + rrii), is positive. Now endow each quark flavor with a new degree 
of freedom, called taste, ranging from 1 to Nt, and replace the Dirac operator D + rrii hj 
[D + mj) ® 1, where 1 is the Nt x A^^ identity matrix.^ Since, evidently, 

det^/^' {{D + nii) ® l) = det(D + m^) , (8) 

where we take the positive Nf^ root, it follows that 



Nf 



Zcont 

/ VA exp{-Sg) Y[ det^/^' (^{D + rrii) ® 1] . 
^ i=i 



(9) 



{rriT^L <C 1), for which a non-perturbative calculation is needed in the effective theory. In that case, 
a non-perturbative generalization of the replica trick is needed to incorporate rooting into the effective 
theory; development of such a generalization is still in progress [xot - 

^ The odd-Nf, continuum negative-mass theory can be reached, at least in principle, by adding an explicit 
theta term with 9 = n to the rooted staggered action. While numerically difficult in QCD due to the sign 
problem, this approach has been shown to work beautifully in the Schwinger model |23j . 

® The "reweighted theories" [7[ provide a way to give this limit a rigorous meaning. The result is a 
formulation with a more complicated gauge sector than in the following discussion, but this difference 
does not affect the fermionic symmetries being considered. See also the discussion of "valence rooting" in 
Ref. 01 . 

^ For the continuum limit of staggered fermions, Nt = 4, but the discussion in the continuum holds more 
generally. 
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Clearly, the combined effect of the two modifications is that nothing at all has changed. The 
Nj:^-Toot theory derived from the A^^t-taste quark fields is nothing but the original theory in 
disguise. Yet, the extra redundancy provided by the new taste index provides a source for 
apparent paradoxes. 

The paradoxes all have their origin in the following fact. The theory with NfNt quarks, 
whose N^^ root is being taken, has, nominally, a much bigger symmetry group than that 
of the original theory. To see this explicitly one can add a (meson) source to the fermion 
determinant [Sj], giving 



Here the massless Dirac operator D and the mass matrix M = diag(mi, m2, . . . ,mNf) carry 
spin, color and flavor indices but no taste index. The source is a color singlet, and has a 
general spin, flavor and taste structure. Because of the taste structure of the source, taking 
the Nj^^ root is not entirely trivial. But since the source is used solely to generate correlation 
functions, we are always allowed to expand as shown. Thus, as before, we only need the Nj^^ 
root of the determinant of the taste- invariant operator {D + M) (g) 1, which is analytic and 
positive. 

The net result, Eq. ffTTj) . shows that one obtains correlation functions of propagators 
evaluated with the correct QCD determinant. One can now derive Ward identities in the 
enlarged fiavor-taste space by doing transformations on J. For example, if all masses are non- 
zero and no two masses are equal, there is a U{Nt) taste symmetry for each flavor, leading 
to NfN^ conserved vector charges. This should be compared to the Nf conserved charges 
in the original theory. If the mass of one particular flavor goes to zero, e.g. — > 0, then 
the up-sector taste symmetry enlarges to an SU{Nt)R x SU{Nt)L x U{l)v non-anomalous 
chiral symmetry, while in the original theory there is still only one conserved charge in the 
up sector. The extra symmetries of the rooted theory mean that there are many extra 
Goldstone bosons (pseudo-Goldstone bosons in the massive case) compared to the original 
theory. This seems paradoxical, since the rooted version, Eq. Q, was supposed to be 
completely equivalent to the original theory, Eq. ([7]). 

However, the extra Goldstone bosons do not modify in any way the physics of the original 
theory. To see this restrict the source to be a taste singlet, J = J 1, where J carries no 
taste indices. Then 



so all correlation functions generated by J are those of QCD. From the left-hand side of this 
equation we see that we can obtain the correct QCD physics in the extended (rooted) theory 
if we restrict ourselves to correlation functions in the taste-singlet sector. Extended charge 
conservation implies that, if the initial state has a total taste-charge equal to zero, then the 
same will be true for the final state. The taste-singlet sector of the extended theory is the 
unitary, physical subspace for QCD with Nf flavors. 




(11) 



(10) 



det^/^' (^[{D + M) (g) 1] + = det^/^' (^[D + M + J] ® 1^ = det (^D + M + , (12) 
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Nevertheless, if we choose to calculate the correlations generated by J in the extended 
theory, using the middle form in Eq. f|T2l) . there will be contributions from taste non-singlet 
(and therefore Goldstone) pions in some diagrams, which perforce must be canceled by 
similar contributions in other diagrams. For example, take J to be a Lorentz scalar and 
flavor singlet, and calculate the contribution of two- meson intermediate states to the two- 
point function. Consider the diagram where the "valence" quark lines (those coupled to 
J) contract in a single quark loop and there is an additional single "sea" quark loop (from 
the determinants). This diagram includes contributions from all {NtN — 1 pions. It is an 
amusing exercise in chiral perturbation theory to check that the contributions of the taste 
non-singlet pions cancel when all diagrams are added. Of course, Eq. (fT2|) tells us that they 
must. 

In the rooted staggered theory at non-zero lattice spacing, the cancellations of the taste 
non-singlet pions will no longer be exact, because the violation of taste symmetry lifts the 
degeneracy of these pions. In particular, there is only one true (pseudo-) Goldstone boson 
per flavor. The lack of cancellation means that the theory is not unitary in the physical, 
taste-singlet subspace at non-zero a. However, it is now clear that the cancellations must 
become exact and unitarity must be restored in the continuum limit. This is not a case of 
"rather subtle cancellations which have not been demonstrated," as Creutz claims. Rather, 
it relies solely on the restoration of taste symmetry in the continuum, for which there is 
extensive evidence, as detailed above. For detailed examples of how unitarity is recovered 



in the continuum limit, see Refs. la, 24 



We note here that the above framework, and in particular Eq. flTT]) . can be used to derive 
non-singlet Ward identities of the type discussed by Creutz in his original posting [2^ (but 
not in Ref. [ll[). Since the taste-singlet sector of the extended, rooted, theory is equivalent 
to the original, target theory, the existence of such Ward identities cannot invalidate the 
rooted staggered approach. 

A question that might have been raised by the discussion so far is the following: Numerical 
simulations of pion and kaon properties with rooted staggered fermions usually use taste 
non-singlet states [H, 0, H, S]- How is this consistent with our statement above that one 
obtains the correct physics with taste- singlet states? The point is that, in the continuum 
limit, the extended vector symmetries show that flavor non-singlet physical states (such as 
the charged pions and the kaons) have identical properties to taste non-singlets with the 
same flavor quantum numbers. So we can choose the taste that is most convenient. This 
is not true, however, for any state which can mix with flavor singlets (such as the ttq when 
mu ^ ma). 

4. Something that Creutz finds "peculiar" is the fact that an ultraviolet limit (a — ^ 0) and an 
infrared limit (m 0) do not commute with staggered fermions. This kind of phenomenon 
is not peculiar to staggered quarks, but is in fact common whenever the lattice theory has less 
symmetry than the continuum theory. For example, in the pure-glue theory on the lattice the 
continuum 2"*"+ glueball splits into two states, the E~^^ and the T^~^, because of the breaking 



of the continuum rotation group to the cubic group. Typically, 112^++ < Tnj.++ [27[, and the 
mass difference is of order a^Ag^^, vanishing in the continuum limit. Because the ultraviolet 
regulator has produced a mass splitting (as it does for staggered fermions), this can lead 
to infrared effects. Define GE{t) and GT2{t) to be (zero-spatial-momentum) correlators of 
interpolating fields in the E~^^ and representations, respectively, normalized so that 
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the lowest states are created with weight 1: 



GE(t) = exp{—mE++t) + ■ ■ ■ ; GT2{t) = exp(— mj,++t) + ■ • • , (13) 

where ■ ■ ■ stands for the contributions of higher-mass states. Then the long-distance (in- 
frared) limit t ^ oo does not commute with a for the ratio of correlators: 

lim lim ^ = 1 ; lim lim ^ = . (14) 

This is the kind of phenomenon that Creutz calls "absurd behavior for a regulator": A 
physical observable (the ratio of correlation functions) develops an infinite derivative with 
respect to the cutoff at a particular point (t = oo). Of course, once one understands what is 
going on, such behavior presents no problem in practice. It is straightforward to extrapolate 
the masses in both channels to the continuum limit to obtain the physical 2+"*" glueball mass. 
Indeed, we never work at the point t = oo on the lattice, just as we never work at the point 
m = 0. 

5. The arguments in this Comment make clear that, in order to use staggered quarks to 
compute physical quantities with controlled systematic errors, it will generally be necessary 
to have a detailed understanding of staggered discretization effects and to extrapolate lat- 
tice data in lattice spacing and quark mass with chiral forms that reflect these discretization 
effects. This is certainly true at lattice spacings of existing ensembles (~ 0.1 fm) where 
the discretization effects are not negligible. The needed forms are provided by staggered 
chiral perturbation theory (SXPT) |22|. With SXPT, even rather complicated physical sit- 
uations such as the scalar correlator (which is sensitive to the anomaly) can be fit and 
properly extrapolated to extract continuum physics [24]. Creutz claims that "completely 
non-perturbative" quantities such as the baryon mass in the chiral limit are "particularly 
suspect." We see no fundamental difference between baryon masses and pseudoscalar de- 
cay constants, which are also completely non-perturbative in the chiral limit and have been 
accurately computed with staggered fermions [3i]. It is true however that staggered compu- 
tations in the baryon channel have been hindered by the absence of SXPT for baryons, which 



has only recently been worked out [26|. With SXPT, well- controlled baryon computations 
with staggered fermions are within reach. Of course, we do not claim that rooted staggered 
fermions are the best approach to calculate all physical quantities in QCD. Depending on 
the size of the discretization errors and how well they are understood, different fermions 
discretizations will produce results for different physical quantities with systematic errors of 
varying sizes. We argue, however, that nothing in Creutz's discussion raises doubts about 
the correctness of properly extrapolated staggered results. 
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